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We determine the behavior of the critical temperature of magnetically mediated p-wave supercon- 
ductivity near a ferromagnetic quantum critical point in three dimensions, distinguishing universal 
and non-universal aspects of the result. We find that the transition temperature is non-zero at the 
critical point, raising the possibility of superconductivity in the ferromagnetic phase. 
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Recent experimental work has shown that supercon- 
ductivity in strongly correlated electron systems is closely 
associated with proximity to magnetic quantum critical 
points suggesting it is mediated by critical spin 

fluctuations [Q, as indicated by theoretical calculations 
HU- However, the interplay of superconductivity and 
criticality is not fully understood. In this paper we study 
the theoretically simplest case, namely p-wave supercon- 
ductivity near a ferromagnetic quantum critical point in 
dimension d — 3. Our work is complimentary to that of 
Abanov, Chubukov and Schmalian j?j who studied pair- 
ing near a two dimensional antiferromagnetic quantum 
critical point. 

We have two motivations. One is practical: one would 
like to know whether the superconducting T c vanishes as 
the magnetic critical point is approached (as shown for 
example in the left-hand panel of Fig. |l] and as found by 
Abanov, Chubukov and Schmalian), or whether it does 
not (as shown in the right-hand panel). The latter sce- 
nario raises the interesting possibility of the coexistence 
of superconductivity and magnetism. This question has 
not been definitively theoretically settled, because nu- 
merical difficulties have prevented a straightforward at- 
tack 0. 
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FIG. 1. Two possible scenarios for the emergence of a su- 
perconducting state near a quantum critical point of a mag- 
netic system: a) the superconducting T c is zero in the quan- 
tum critical point (QCP), and alternatively b) T c is finite in 
QCP. The parameter r measures the distance from the quan- 
tum critical point r c [ pd|Jl^ ]. In experimental realizations r 
corresponds to hydrostatic pressure 

Our second motivation is theoretical. Studies of mag- 
netically mediated superconductivity have almost uni- 
formly been based on the Eliashberg equations (defined 



below) Q which are simple generalizations of the equa- 
tions which describe conventional phonon-mediated su- 
perconductors ||. While these equations are believed 
|l0f to give the leading contributions to the low-energy 
behavior of systems near critical points, non-critical and 
high-frequency processes may also be important for the 
superconducting T c (These lead, e.g. to the pL* famil- 
iar from conventional superconductivity) |^| . Our results 
show how to isolate the critical contributions and allow 
the magnetic analogue of fi* to be estimated. 
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FIG. 2. Eliashberg equations in diagrammatic form. 
E(p, iuj) — iu{l — Z p (u)) is the normal-state quasiparticle 
self energy and W(p,iui) is the anomalous self energy. 

We consider a three-dimensional metal near a ferro- 
magnetic quantum critical point. The magnetic suscep- 
tibility is (n]JT| 
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where r is a parameter that measures the distance of the 
system from its quantum critical point (Fig. [l]) and the 
ellipsis denotes less singular terms. Here pf is a mo- 
mentum scale of the order of the Fermi momentum and 
A ~ vfPf is an energy scale of the order of the Fermi en- 
ergy. We assume (following ||||,|],[lO|] ) that the coupling 
of these to the electron system is given by the Eliash- 
berg equations (shown diagrammatically in Fig. ^) for 
the electron self-energy £(p, iui) = iu){l — Z p (uj)) and the 
anomalous self-energy W(p, iuo). W(p, iuj) vanishes at the 
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superconducting critical temperature and grows continu- 
ously in the superconducting state. We find T c by solving 
the linearized Eliashberg equations, which are 

Z p {lo)) = g 2 s [ N{n p ,)AVL p , f de p > (2) 

J J — oo 

^tV] X (p-p',iuJ-iu}') , 1 T 

w(p,iLo) = g 2 si { N{n p ,)&n p , { dt p , (3) 

J J — oo 

x (p-p,%u>-%w) - 2 -. 

[tu'Z p {u')] -t z p , 

Here the momentum integration has been separated into 
integration in a direction perpendicular to the Fermi sur- 
face (e p integration) and integration over angles £l p of the 
spherical Fermi surface; N(fi p ) is the angle-dependent 
density of states of the quasiparticles on the Fermi sur- 
face. The numerical factors si relate to the nature of the 
spin fluctuations and the symmetry of the pairing state. 
For a system with a Heisenberg symmetry there are three 
independent soft spin components all of which contribute 
to loZ so Sq = 3. However, for spin triplet pairing only 
one combination can contribute to any given component 
of the gap function, so si = 1. The importance of this 
factor was stressed by Monthoux and Lonzarich For 
a system with a strong Ising anisotropy, both sq and 
,si = 1. We will present results for the Heisenberg - Ising 
crossover elsewhere, g is a constant vertex representing 
the interaction between spin fluctuations and low-energy 
quasiparticles. It may be experimentally defined from 
the singular (as r — > 0) behavior of the specific heat co- 
efficient 
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Eqs. (g) and (g) apply only for frequencies much less 
than the electron bandwidth and only if the momen- 
tum dependence of Z p and W is negligible relative to 
the frequency dependence, conditions which are satisfied 
for the leading singular behavior as r — » 0. We therefore 
employ the Migdal approximation Jl3[ Z p {ui) — ► Z(u>), 
W(p,oj) — ► W(ttp,iui) and perform the integral over the 
magnitude of the momentum. To perform the remaining 
integration over angles we note that iu>Z(u>) has the full 
symmetry of the lattice, while for p-wave superconduc- 
tivity W corresponds to the / = 1 spherical harmonic. 

The momentum transfer q carried by the spin fluc- 
tuations in Eq. (Q) is given by q 2 = (p — p') 2 = 

2p% (1 - (p ■ p')/\p\\p'\) + 4'/ v f- The first term in I 2 
is obtained by placing both momenta p and p' on the 
Fermi surface while the last term is a small correction 5p 
taking into account the fact that intermediate states can 
explore regions close to the Fermi surface (Fig. ||) and 



will be important as a cutoff. We perform the e p integral, 
use the angle dependences of Z and W and obtain 

M (1 - Z(u)) = -nTs J2 D a(u- w')Bgn(w') (5) 
W, M = nT Sl Y,Di{u-w>) Wl }f ] m , (6) 



where 
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with U = [(v/x)t&iT 1 (x/v) + x 2 +r] , v = u/A, 
Pi{x) is a Legendre polynomial and the \u>Z\/A comes 
from the e 2 , term. It is numerically very small and is 
important only as a cutoff at r < T 2 and v = 0; except 
in these cases we drop it. In the following we combine all 
constant prefactors in Di into a single coupling constant: 
l6ir 2 g 2 N -► A. 



|5p|= Ep ./v F 




FIG. 3. Fermi surface with momenta participating in the 
interaction. The dashed line is the momentum transfer q when 
the correction Sp is neglected. 

To solve Eqs. (@) - @ we follow Bergmann and 
Rainer JTJJ, defining a new order parameter $;(w) = 
Wi(iuj)/\ujZ(Lu)\ and casting Eqs. (||) and (||) into an 
eigenvalue problem for an eigenvalue p 



E 



|q/Z(a/)| f 
siDi{u - u)')- s Q — 'Su. 
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$j(</) = p$j(w) (8) 



\luZ(lu)\ = \lu\+ttT I D (0) + D o{" ~ £■/) ) • (9) 

At high temperatures the eigenvalues p n (T) are nega- 
tive; at T c the leading eigenvalue crosses 0. We solve the 
matrix system numerically; the size of the kernel, K nm — 
siDi(uj n - uj m ) - S nm s \oj m Z(u! m )\/irT is - A/(2irT c ). 

For p-wave pairing in systems of Heisenberg symme- 
try the critical temperatures are typically irT c ~ 10 _5 A 
which translates into numerically unmanageable kernel 
sizes of N ~ 50 000. We therefore use a down-folding 
procedure: we separate $z(w„) in Eq. (N) into a low- 
frequency part &f ow (uj n ) with < |n| < Nlow and 
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high-frequency part &f IGH (ui n ) with Nlow < \n\ < 
N. Then Eq. (||) can be written as a block lin- 
ear system and formally solved for $> HIGH ; yielding 

K LOW . $LOW = p$ IOW with J(LOW = + 

Z)|i|,|j|>iViow "^ n< ^ _ ^')~ ^3»n- This transforma- 
tion is exact. The simplification is that for large Nlow 
K is nearly diagonal so K^ 1 may easily be computed 
in the "high" subspace. In physical terms, this approx- 
imation retains only the diagonal scattering-dominated 
terms K nn ~ {s x - s o )D o (0) - (2n + 1)(1 + (2/3) Xs + 
(2/3)Aso ln(A/27rTn)) and drops all off-diagonal pair- 
ing terms siD\(v nm ) ~ (Asi/3) ln(A/27rT|n — m|) for 
n, m > Nlow m the high-frequency kernel. We have 
verified that this approximation reproduces faithfully the 
eigenvalues of Eq. (||) for large temperatures, and that 
our results are insensitive to the choice of Nlow- 

Our results for T c (r) are shown in Fig. ||. The top 
panel of the figure demonstrates the convergence of the 
scaling procedure with reduced kernel size Nlow- Ker- 
nel sizes Nlow > 500 show satisfactory convergence, so 
we have used sizes Nlow = 500 in most of our work. 
That large kernels are needed shows that in this problem 
T c is not controlled by low-energy physics. As previously 
noted (|] T c is very low in the p-wave Heisenberg case 
because of the factor of three between the pairing vertex 
and the self energy. The Ising case has not been previ- 
ously studied; we see T c is much higher. 

We find that in both Heisenberg and Ising cases, 
T c (r — > 0) > 0, raising the interesting possibility of su- 
perconductivity extending into the magnetic phase. We 
confirm that T c (r = 0) > using a variational argu- 



ment. The ansatz W{u> m ) — AZ(w m )0(w* 



J al- 



lows the leading eigenvalue to be computed if to* « A. 
At r = 0, the leading eigenvalue becomes positive for 
T < T var with 



T var = _ exp 
7T 



3/2A + s so 
' si (1 + ln(A/w*)) In 



(10) 



which is thus a lower bound for T c . Abanov, Chubukov 
and Schmalian jTj , who studied a two dimensional antifer- 
romagnetic problem, argued that the divergent mass en- 
hancement associated with the critical fluctuation would 
drive the superfluid stiffness and thus T c to zero. In their 
case the divergent mass occurs only at one point on the 
Fermi surface, so it seems to us the considerations of Hlu- 
bina and Rice []l5| should imply a non-zero superfluid 
stiffness. In any event, in the ferromagnetic problem of 
interest here the critical fluctuations are long- wavelength, 
and thus do not lead to divergences in the "transport 
mass" controlling the superfluid stiffness. 
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FIG. 4. Results for T c (r) a) near r — for a set of ker- 
nel sizes Nlow and b) in a broader range for two coupling 
constants A = 1.5 and A = 10.0; in the Ising case (so = 1) 
we have plotted T c /100 for better visual comparison with the 
Heisenberg (so = 3) curves. A ~ 2pfVf is the characteristic 



spin-fluctuation frequency; the curve 2-kT c — 
two regimes when T c < r 3 ^ 2 and when T c > 



r 3/2 
„3/2 



separates 
For the 



Heisenberg case (so = 3) Nlow 
case (s = 1) N = A/2ttT c . 



500 while for the Ising 



We see from Fig. [|a that the Heisenberg case T c (r) 
displays a maximum at a small non-zero r, whereas in 
the Ising case there is no maximum. We believe the small 
r behavior is controlled by the interplay of pairing and 
scattering, as discussed by Bergmann and Rainer [ fl4| for 
s-wave superconductivity and by Millis, Sachdev, Varma 
for d-wave 01|. To see this mathematically we compute 
dT c /dr using the Feynman-Helhnan theorem 113]: 



dr 



dp_ 



'dp 
dr 
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dr 



(11) 
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The expectation value is infrared dominated numerically 
and we find that at small ui, $>i(uj) ~ From Eqs. 

(0) and (||) we see there are two contributions to dK/dr: 
one positive and proportional to si , coming from the pair- 
ing term Di(y) and one negative and proportional to so, 
from the depairing term \u>Z\ in Eq. (|8|). As r — > the 
dominant term in Do becomes identical to the dominant 
term in Dy. It is convenient to isolate the contribution 
from zero- frequency spin fluctuations. For the leading 
singular behavior in r we find 



by the low energy spin fluctuations which are theoreti- 
cally tractable, and demonstrated the crucial role played 
by the symmetry of the magnetic fluctuations. We have 
found generically that T c > at the magnetic critical 
point, raising the interesting possibility of superconduc- 
tivity within the ordered phase. 

We acknowledge NSF DMR 9996282. A.J.M. thanks 
G.G. Lonzarich for many helpful discussions. 
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2Vr(Vr+KZ|/A) 
comes from differentiating Eq. (0) at i> = while 
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(14) 



comes from differentiating Eq. ([?]) at v ^ and dropping 
the |u;Z| term. The scaling function 



F(x) = / 
Jo 



ydy 



(z/y + y 2 + l) 2 



(15) 



F(0) = 1/2 and as x -> oo, F(x) -> (27r/9V3)a;- 2/3 . 

For sq = 1 (Ising case) the 2?o term vanishes and the 
term is negative. The pairing and depairing effects 
of quasistatic (u> < T) spin fluctuations exactly cancel 
(as in the s-wave case Q) while at uj > T the pairing 
effect wins. Thus T c monotonically increases as r — > 
because the spin fluctuations become stronger. At r = 0, 

dT c /dr T- 2 / 3 , i.e. T c {r) is linear; for r > T c 3/2 the 

derivative — > (lnl/r)/r, so we expect T c ~ In 2 1/r. 

For so — 3 (Heisenberg case) the T>q term is non- 
vanishing, and indeed is dominant at small r: quasistatic 
spin fluctuations are pairbreaking. At r — T c is set 
by the temperature at which the effect of these fluc- 
tuations becomes small enough to allow pairing. For 
r < \ 2 {ttT c ) 2 In 2 A/T c , the luZ term is important and 
dT c /dr - l/(V? T A7rT c lnA/r c ). For \ 2 (irT c ) 2 In 2 A/T c < 
r < (7tT c ) 3 / 2 ; dT c /dr - 1/r. In our calculations, this in- 
termediate regime is not wide enough to see. For larger 
r, the variation of the pairing (2\,) term with r becomes 
most important. 

To summarize, we have presented a theory of the vari- 
ation of a p-wave superconducting T c near a ferromag- 
netic quantum critical point. We have shown that the 
variation of T c with distance from criticality is controlled 
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